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The approach supposed by Feynman [l| for disentangling of non-commuting operators (FDM) was successfully 
applied to various time-dependent problems in quantum mechanics _ @ ■ Using FDM one can [J 0| calculate the 
transition probabilities between eigenstates \m,t — > — oo) and \n,t —* +00) of a quantum oscillator with definite 
initial and final quantum numbers, and compute corresponding generating functions. If the oscillator has a constant 
frequency lu being under an external force /(t) , for the generating function one has 2 : 

G{u,v\u)= V w mn {u)u m v n = (1 - uv)- 1 exp \-u (1 ~ v) } , (1) 

* — ' "I 1 — uv 

m,n—0 v ' 

where u, v are the auxiliary complex numbers and v is the oscillator excitation parameter: 

1 



V = 7T- 

2ui 



f{t) exp(—iu>t)dt 



, h = m = l. (2) 



From ([T]) one obtains some non-trivial relations for probabilities w mn . First, integrating by v, one can derive several 
peculiar sum rules 



W m n(y)dv = 1, (u) mn = J W mn (v)vd,V = TU + U + 1, 

(3) 

(Af 2 ) mn = / w mn (v) [v - {v) mn \ 2 dv = 2mn + m + n + l 



for arbitrary quantum numbers m and n. For the initial vacuum state (m = 0) we have (Av' 2 }o n = (f)on which 
corresponds to the Poisson distribution. 
Taking in |T]) u = v one has 

G(u, u\v) = (1 - u 2 )" 1 exp {2vu{\ + uy 1 - v) 

which yields 

S k{^)= ^2 w ™-n{v) = e~ U Pk(v), (4) 

where pk{v) = X)s=o( — ^YL s {2v) and L s are the Laguerre polynomials, so 

Po = l, Pl = 2v, p 2 = 2v 2 ~ 2v + 1, Pi = ^ 3 - Av 2 + Av, ... 
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1 A natural condition f(t) — » at t — » ±oo is supposed. 

2 This expression can be easily derived 0] from Schwinger result [f| for a sum X]m=o w mn u 
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For the oscillator with variable frequency u>(t) without an external force 3 (f(t) = 0) one obtains 

oo 

G(u,v\p)= w mn (p)u m v n = v/(l - p)j [(1 - uvf - p(u - v)% (5) 



m,n— 

and 

l 



2 

G(u,v I p)(l — p) dp = (Arctanhu — Arctanhw). (6) 

u — v 



o 



Here p, < p < 1 is the oscillator excitation parameter; details can be found in [l, Q. The frequency uj(t) is an 
arbitrary real time function. As usual, we propose the boundary conditions: 

uj(t) — > ui± at t — > ±oo 

which allows one to define the final and initial eigenstates of the oscillator. Note that the expressions ([1]) and ([5]) 
result from the general Husimi expression 0]. 
As a result we have 

1 w mn { P ), i + (-ir+" 

dp — , 

1 — p m + n + 1 

(7) 



o 

l 



«Up), 1 + (-!)"+" 

— 7r= c ^ = — i i - ' m ^ n - 

/VI _ P | m-n | 



i 

Analogously to ([3]) let us calculate the integral J m „ = J w mn {p)dp. For diagonal (m = n) transitions 

o 



Jn 



2n+ 1 



1 + 



(2n + 3)(2n- 1) 



, n = 0,1,2,..., (8) 



at to n the expression for J mn is more cumbersome. 
Finally, taking in ([5]) u — v, one obtains 



G{u,u\p)- v ' 



1 — it 
and 



K/T— 5, fc = 0,2,4,... 

^•(p) = w "^p) = \ ( 9 ) 

m+n=k [^0, fc= 1,3,5,... 

(compare to (|4])). 

On differentiating subsequently the generating functions on parameters u and v, one can calculate the average 
quantum number in the final state 

dG / \ m / \ 1 , , 1 \ 1 + Z 9 

0=1 m=0 n=0 

the dispersion (An 2 ) m and other higher distribution momenta 



Y] (n) m u m , (n) m = y nw mn = -- + (to + - ) , (10) 



For this case the transitions occur only between states with same parity. 
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In [4| a more general model of a singular oscillator with variable frequency was considered: 

zj 1 2 i 1 / j. \2 2 i .9 

ff=- P + - W (t)* + (n) 
< x < +oo, g — const, g > — 1. 
It is well known that at a fixed t the instantaneous spectrum of the Hamiltonian (fTTj) is equidistant (see, e.g., (§1): 

E n = 2u(n-j), J = -\-\^/T+~g, n = 0, 1, 2, . . . . (12) 

Using FDM, one can disentangle the operators in H and, after application of some group theory methods, calculate 
the transition amplitudes between initial \m) and final \n) states, being expressed in terms of the generalized Wigner 
function for the irreducible representation of the su(l, 1) algebra with weight j. As a result, one obtains explicit 
analytical expressions for the transition probabilities w mn . The corresponding generating function takes the form: 

g{u, v) = V w mn u m v n = A 2 \ 2 (13) 
' 1 — uu ■ A^ 

m.n 

where 

2(1-/0 



1 — p(tt + u) + to + \J [1 — p(it + v) + uv] 2 — Auv(l — p) 2 



(see @ for details). 

One can show that at j = —3/4, which corresponds to the regular oscillator (g = 0) with odd levels, these expressions 
turn into 

The most interesting physical case is the excitation of the oscillator ground level, i.e. m = 0. Taking in (|13[) u = 
and, correspondingly, A = (1 — p)/(l — pu), one obtains 



sM = £« n =(^J (w) 



-2j 



and 



^o»(p) = ^ r ^p"(l-p)- 2i , n = 0,1,2,... (15) 

This formula also yields Q the transition probabilities of the regular oscillator for j = —1/4 (even states, th level 
excitation) and for j = —3/4 (odd states, 1 th level excitation). 

To conclude, let us note that the generating function (| 1 3[) determines also the adiabatic expansion of the probability 
w m n for the case of small values of the excitation parameter. For example, for diagonal transitions (to = n) one has 

w nn {p) = 1 - 2[n 2 - (2n + l)j]p + ■ ■ ■ = 1 - \ (N 2 + N + l) p + 0(p 2 ) (16) 

where 



N = 2V(n - j) 2 - U(J + 1) + 3/16) - -. 
In particular, N — 2n for j = —1/4 and N = 2n + 1 for j = —3/4. 
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